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CONCAVE  .►^ncRmnS-  KCS  GAS'.LJ  C  3LE^'DS 


A.  S.  *ianx>e 


I ,  Tnt ro-1  jctlon 

-.I-  ;  ''pc-r  presents  a  -oJel  cf  the  ecanoric?  of  notor  (rascline  tlcr.l- 
Ir,'  —  a  T.^xt.-2l:lnc  probltr  suijtct  to  certain  r»  strain*;  .  -  — 

of  tae  '••'ill  c«;».';ot  te  i.t-.llcj  <.t:.in  t*i'»  IraTf-*--'  f"  ^ 

Lafirangean  multlpllere.  As  yet,  nelth“r  George  Dantzlg  nor  I  nave  dj:- 
ceeied  In  converting  the  tysten  into  a  linear  prOtra-T.-iinr;  ret'v 
is  short  enough  to  be  convenient  for  coeputa ticnal  purposes,  ml  A.’.icn 
at  the  saae  time  preserves  tne  octane  technology-. 

Kuhn  ni  Tucker  hare  e.-’tablished  the  u.eorei  that  a  scluticr. 

to  *Jnig  type  of  maximizing  problem  Is  equivalent  to  a  certain  sinimajt 
sadilepoint,  and  at  the  sjggcstivn  cf  Harry*  Vark  'it:,  the  rasolme 
blenling  model  was  recast  into  tne  appropriate  minimax  form.  An  iter- 

y 

atire  dirital  tecimique  war  tn*n  -ievlcei  I'jr  solvlnr  the  mlrjimax  ir-  1*5T, 
From  a  rigoro’JS  matne-itica]  stanJpoi-it,  tne  f  resent  Iteratl-^e 
technique  ir  by  no  means  s.i‘i  f.i.tcry.  The  net.^i  ;.r«s'.tr,  lead 

tc  S''iutl''n3  that  rhu];  t «’  sif  i"!  clcntly  accurat'*  fren  *hf.  r'^andpolnt 
of  refinery  -nanarfrF.  lonvi  rf"  •..c  taker  place  witt.out  an  Inorilna^e 
amount  of  cnmp^.tirj’  tim*  .  All  cilc'jiatior.s  eire  pvrl'-imni  on  a  r.a:  .  :i( 
that  is  realily  available  tiir^ur  out  the  '  nittd  Gistc'  --  an  liV  Card 

1/  I  should  li’^e  to  ncknowlc dj-e  ny  indel  l»-dne5e  to  Harry  Uarkceit.:, 

Ge  :  ye  Dariitiy,  and  to  ItnviJ  Lmpfield  for  their  many  contriboti.  no 
'  t.-ti;  tlroc  of  rtoearch. 

r/  Varkowlt:  is  cjrr*r.tly  oxperlmf  ntirg  wl  tli  an  analcgue  r.achl.ne 
”  F. IrLuax  solution. 


'..etu,'’ 


PrDgTAaaed  Calculator,  Model  ?.  A  general  pjrpcse  floatlnf.-  :i»  .  :<il 
was  employed.  On  the  first  problem  run  off,  approximate!;/  te:.  r.o'ors  of 
running  time  were  required.  For  the  sixth  problem  cf  this  type,  the  running 
time  had  been  cut  down  to  less  than  three  hours. 

In  order  to  maintain  a  check  on  the  solutions,  the  mo-del  war  kept* simple 
encjfh  30  that  it  could  also  be  solve!  by  gra:.ntcal  techniques.  For  that 
reason,  the  model  will  not  in  Itself  cc  of  interest  to  refinerv'  op-rators. 

From  tholr  standpoint,  nevertheless ,  tne  computatlo:is  art  sipnlflcarit.  The 
results  bring  into  question  the  validity  of  a  number  of  rulcs-of-thjmb  that 
are  current  In  the  industry  for  solvlnr  this  type  vf  econcriic  problem.  A 
second  point!  the  calculations  five  s  -me  indication  of  the  financial  licprove- 
r.ent  that  may  be  obtained  by  takinK  account  of  non-llnoar  octane  blenllng 
relationships.  And  finally,  these  machine  runs  nave  serred  as  a  pilot  study, 
prellminajy  to  attacking  problems  that  Involve  larger  numbers  of  variables, 

Tr.e  material  will  be  presented  in  the  following  sequence:  octane  miaber 
prediction  devices;  the  economic  model;  the  computational  procedxire;  and  then 
a  Jiscusrion  of  tne  results  of  the  six  CPC  runs. 

?.  Octane  number  predictions 

Like  most  kinds  of  engineering  work,  the  forecast  of  gasoline  octane 
numbers  is  more  cf  an  art  than  an  exact  science.  Refiners  themselves  neces¬ 
sarily  make  paper  forecasts  of  the  octane  number  of  proposed  blends,  but  before 
marketing  a  product  they  will  almost  inevitably  take  the  precaution  of  testing 
ti.e  mixture  In  in  octane  ratinr  engine. 

(a)  Tetraethyl  leal  concentration  levels. 

One  phase  nf  the  problm  nas  been  Investigate!  extensively  by  the  refining 
l.ndustry  —  the  relallonehlp  between  octane  nunoerc  and  the  amount  of  tetraethyl 
lead  added  por  gallon  of  a  particular  grsollne.  It  hae  been  invarlahly  ■'bsepveJ 
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Fig.  I  — Octane  number  versus  TEL  concentrotion  for  three  gasolines 
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lh*t,  as  the  lead  concentration  level  Is  Increased,  the  octane  number 

increases  at  a  decreasing  rate.  In  order  to  predict  this  relationship, 

there  Is  In  widespread  use  an  ethyl  blending  nomogram,  reproduced  as 

Figure  The  points  obtained  by  testing  a  particular  gasoline,  when 

plotted  on  such  charts,  will  generally  lie  along  a  straight  line.  The 

line  labelled  "1"  represents  three  observations  for  a  particular  gasoline 

1/ 

reported  by  3ogen  and  Nichols.”  The  lines  labelled  "2"  and  ”3"  represent 
two  more  sets  of  observations  reported  by  them.  In  order  to  cut  down 
on  the  nuaber  of  knock  engine  runs,  It  Is  the  general  practice  to  test 
one  sample  of  clear  gasoline  (l.e.,  no  TEL)  and  another  sample  of  the 
same  gasoline  containing  3  ml.  of  the  ethyl  fluid.  (3  ml.  typically 
represents  the  maximum  allowable  concentration  in  American  motor  gaso¬ 
lines.)  These  tests  yield  two  points  which  can  be  plotted  on  the  blending 
paper.  The  octane  numbers  for  intenaedlste  lead  concentration  levels 
can  then  be  predicted  by  straight  line  interpolation. 

The  eti^l  blending  chaa*t  that  is  In  general  use  —  the  one  put 
forward  by  Hebl,  Rendcl,  and  Oarton  In  1939  ^  ^  —  was  itself  derived 
frean  empirical  observations,  and  not  from  any  algebral?  formula.  I 
have  found,  however,  that  there  Is  a  three-parameter  analytic  expression 
that  gives  a  close  approximation  to  the  results  predicted  by  the  blerding 
chart.  I'lhere  t  represents  the  octane  number  of  the  leaded  gasoline,  and 
X  represents  the  TEL  concentration  level  (in  ml.  per  gallon) i 

(1)  t«a»bxe  c 

1  ♦  X 

l-cr  unki.own  gasoili>e,  the  tiirec  parameters,  a,  b,  and  c,  have  tc 
t  i tier-'?. red.  This  may  be  accomplished  rtallly  by  taking  ihree 


r--3o3 

“5“ 

observations  of  octane  nmber  for  v<irioj*  lead  levels^  and  5olvin^  t'-.rtt  slr.ul- 
♦aneoM*  linear  equations.  Two  observations  can  be  the  expcrinental  rncB»  and 
the  third  =ay  be  deterrined  in  tre  standard  way  on  tne  ethyl  blending;  chart. 

U5lr.g  equation  (1)  in  this  way,  approxiaately  20  seta  of  constants  have  been 
detencined  for  actual  gasolines.  Testing  seven  lead  levels  in  each  of  these  20 
cases,  the  equation  has  virtually  always  given  a  prediction  that  lias  within  * 

0,2  cetane  nutrber  of  the  one  yielded  by  the  blending  chart. 

No  PLCclal  significance  jhoulc  be  attached  to  the  particular  lore;  cf  equation 
(1).  It  is  T.frfly  a  device  for  enabling  a  conputing  •nachlne  tc  perforr.  the  s^-ae 
calculation  as  a  refiner  with  his  blending  chart  and  his  straight  edge. 

(b^  Two-cooponent  gascline  blends. 

A  second  phase  of  the  blending  probleoi  is  more  controversial  than  the  TbL 
asnect.  Kor  a  lead  level  of,  say,  3  nl.,  and  for  a  50-^^C  irdxture  of  two  gas- 
cllnes  A  and  B,  refiners  frequently*  calculate  the  octane  number  of  the  blend 
to  be  the  50-1^0  weighted  average  of  the  octane  number  of  gasolines  A  and  3  ~ 
each  with  3  al.  T.he  weighted  average  appears  ratlsfactory  for  itany  gasolines  -- 
especially  when  both  oorp'onents  are  of  paraffinic  nature.  There  arc,  hoi»ever, 
at  least  tac  papers  publicly  available  —  one  hy  Eastman  7  and  the  other 
by  Bcgen  and  Nichols  "•  that  call  into  question  the  strai  gjit-line  averag¬ 

ing  xethed.  Both  pap-er^  indl.uite  that  as  the  p>crcentagc  of  the  high  octane 
cQupcnent  in  a  binary  nixtvre  Increases  (both  coBponents  initially  at  the  sane 
TEL  concentration  level),  the  octane  level  of  the  mixture  nay  Increase  at  a 
decre asi r.g  rate.  irj cally ,  the  cetane  number  is  a  concave,  nonotone- 

Increaslng  f’uncMcn  of  the  pcrcentarc  cf  the  high-octane  component  in  the 
blend.  In  ordinary  Icniiiage,  the  octa,ne  n\B.Ler  of  the  blend  tends  to  exceed 
the  welrlited  average  of  the  octane  nur.bcr  of  the  two  components. 


-  V. 


7 


•*. 


j*  rjt  ^  r*^ 


•*. 


o 

«>  ee 


”5  06, 


■f  3.0  ml  of  TEL 


♦  I  0  ml  of  TEL 


* 

I  02 

u 

O 


vCleor 


0  20  40  60  00  100 

%  thermolly  cracked  gasoline  in  binory  blend 

Fig.  2— Octane  number  versus  percent  of  thermally 
crocked  component  in  binary  blend.  Three  TEL  levels 
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3o,-en  and  Nichols  are  primarily  concerned  with  three-component  mlr- 
tures,  and  they  report  only  five  observations  at  each  of  three  TEL  levels 
for  strictly  binary  mixtures.  The  set  of  15  points  has  been  plotted, 
and  is  reproduced  as  Flpure  2,  These  15  observations  are  the  basis  for 
constructing  the  present  model.  For  each  TEL  concentration  level  shown 
on  this  chart,  straight-line  interpolation  has  been  used  between  adjacent 
observations.  In  all  likelihood,  these  straight  lines  also  underestimate 
the  octane  nvinber  that  is  attainable  by  blending.  The  existing  data, 
unfortunately,  do  not  provide  a  basis  for  curvilinear  interpolation. 

At  first  glance,  it  would  appear  to  be  a  trivial  matter  whether 
octane  ninbers  were  calculated  by  drawing  strai^t  lines  between  two 
adjacent  points  or  between  two  end  points  of  a  constant-TEL  line.  Never¬ 
theless,  a  sample  calculation  presented  below  in  Section  (6)  indicates 
that  the  two  methods  lead  to  TEL  input  requirements  that  differ  by  20il 

3.  An  economic  model  of  the  blending  problem 

It  is  assumed  that  the  refiner  has  available  fixed  quantities,  and 
q2,  of  two  blendings  stocks  -  of  types  1  and  2,  respectively.  The  ^1 
stock  Is  the  hlgfi  octane,  cataly tlc-cracked-plus-polymer  gasoline,  tested 
by  5ogen  and  fhchols.  The  ^2  stock  is  their  low  octane,  thermally  cracked 
gasoline.  These  two  may  be  blended  together  In  any  proportions,  and  may 
be  mixed  with  tetraethyl  lead.  Ther<'  are  two  products  —  1  and  2  — 
premium  and  regular  grade  gasoline,  respectively.  These  two  are  to  mtet 
minimum  F-1  (Research)  octane  number  specifications  of  and  Nj,  and 
they  may  be  sold  at  unit  prices,  and  p^.  The  cost  of  1  ml.  of  tetra¬ 
ethyl  lead  is  represented  by  p^. 

Instead  of  using  the  two  blending  stocks  directly  foi  producing 
finished  gasoline,  the  refiner  may  also  elect  to  produce  a  fictitious 
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intennedlatc  product  -  "stock  3"  -  by  mixing  together  tf.e  two  initial 
materials  in  70-30  proportions.  This  third  stock  may  in  turn  be  bltria^d 
into  either  of  the  two  finished  products.  By  permitting  the  refiner  to 
prO'Juce  this  fictitious  intermediate  product,  the  model  in  effect  permits 
him  to  interpolate  between  all  adjacent  observations  on  Chart  2.  There  are 
eight  independent  variables  in  the  problem,  defined  as  followsi 


(2) 


m  number  of  gallons  of  stock  1  in  product  1 

r  number  of  gallons  of  stock  2  in  product  1 

x^j^  -  nuaber  of  gallons  of  stock  3  In  product  1 

xjj]^  M  -.1.  of  tetraethyl  lead  per  gallon  of  product  1 

x-^2  *  number  of  gallons  of  stock  1  in  product  2 

X22  z  number  of  gallons  of  stock  2  in  product  2 

X32  ■  number  of  gallons  of  stock  3  in  product  2 

X|j2  =  ml.  of  tetraethyl  leal  per  gallon  of  product  2 


The  set  of  8  xj  will  be  indicated  by  the  vector  x.  The  refinery's  net 
gasoline  realization  will  be  temed  g(x).  The  net  realization  on  one  gallon 
of  product  equals  the  selling  price  of  that  product,  less  the  TEL  cost  per 
gallon. 


(3)  g(x)  r  (pj^  -  (xii  ♦  X21  *■  xji) 

♦  (pp  -  P3Xh2)  (xi2  ♦  X22  ♦  xip) 

The  exprer.slon  for  g(x)  is  to  be  maximized,  subject  to  certain  restrictions, 
Indicated  below  by  (li)  -  (7),  (9),  and  (10).  The  first  condition  is  that 
the  Xj^  must  not  be  negativei 

(I4)  *1  >  0  ff’*'  all  1 
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.•(ext,  the  TEL  concentration  levels  nust  not  exceed  3  ml.  per  fc;allon. 

(5)  3  -  >  0  J  -  1,  2 

The  blending  schedulei  l.e.  the  vector  x*  must  not  call  for  more  of 

the  blending  stocks  than  the  quantities  of  them  that  are  available.  The 

excess  quantity,  if  s.ny,  of  blending  stock  1  will  be  denoted  by  fj^,  and 

the  excess  of  blending  stock  2  by  £2* 

(6)  f^Cx)  =  qj  -  xj^j^  -  AJ2  ■  0*7(x33^  ♦  x-,-)  ^  0 

(7)  f$fx)  •  qj  -  X21  -  X22  -  0.3(x3j^  ♦  X32)  >  0 

Finally,  the  teo  gasoline  products  must  meet  the  appropriate  octane 
mxnber  specifications.  To  simplify  the  notation,  a  new  variable  is 
defined,  t^j.  This  variable  indicates  the  number  of  octane  points  by 
which  stock  1  exceeds  (algebraically)  the  specifications  for  product 
when  xjjj  al.  of  lead  arc  aided  to  stock  i.  The  variable  t^j  Is  evidently 
a  function  of  x^^j,  and  Is  determined  through  the  following  relatloni 

(8)  *  aj  ♦  bi  x,^j  e  _!l -  - 

1  V  X|^j 

In  order  to  ensure  meeting  octane  number  specifications,  the  follow¬ 
ing  conditions  must  then  applyi 

(9)  f3  w  *il*^ll^*i4l^  *  *?l*^21^*ll)  ♦  *31*^31^*U^  >  0 

(10)  fjj  .  xi2.ti2(x|p)  ♦  X22>^?t(^L?)  *  ’'3?*<32f*L2)  t  ^ 

The  mathematical  problem  posed  here  is  the  raaximlration  of  expression 
(3),  subject  to  conditions  (L)-(7),(9),  and  (10).  Before  going  on  to  the 
computational  procedun  ,  it  is  worthwhile  to  reexamine  some  of  the  assump¬ 
tions  that  have  been  slipped  lnv.o  the  analysis. 


A  .A., A  i«.A  ._A  A.A  -A  A.\  r.A  A.A »,  L 
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Alth  the  excep  tion  of  the  octane  blending  rcl at lonsr. ,  t;  eoe  arr/oa,- 
tlcns  closely  restinble  the  ones  underlying  the  avgas  linear  prcfiXacjiing 
model  of  Cnarne*,  Cooper,  and  Mellon  /~2_7.  Like  the  earlier  model,  I  hare 
arcmcd  that  the  quantities  of  blending  stocks  available  and  that  the  octane 
characteristics  of  these  blending  stocks  are  Independent  of  decisions  In 
the  blending  department.  In  fact,  from  the  viewpoint  of  a  refinery  super¬ 
intendent,  these  quantities  and  characteristics  are- variables  subject  to 
control.  By  altering  the  reactor  temperatures,  the  recycle  ratios,  and 
the  asslgment  of  Intermediate  distillate  oil  streams,  the  central  management 
is  able  to  Influence  the  site  and  composition  of  the  gasoline  blending  materials 
This  type  of  interlock  can  only  be  studied  In  a  larger  model  that  cuts  across 
departmental  lines. 

Again,  like  the  linear  programming  system,  the  present  one  assimcs  th.at 
unlimited  quantities  of  tetraethyl  lead  may  be  purchased  at  the  going  rate, 
and  that  unlimited  quantities  of  premium  and  of  regular  grade  gasoline  may 
be  marketed  at  the  current  spot  price.  Such  assuoptlonr.  are  probably  real¬ 
istic  for  the  smaller  refiners,  but  they  are  open  to  serious  doubt  In  the 
case  of  the  major  ones.  For  the  latter,  a  more  complex  type  of  g(x)  payoff 
function  Is  certainly  In  order.  Lacking  this  more  complex  function,  I  have 
fallen  back  on  the  simpler  one.  The  econcolst  can  at  least  comfort  ninself 
with  the  thought  that  equation  (3)  Is  the  one  usually  employed  by  Individual 
dep.-u-tments  within  a  large  refining  organization. 

In  one  respect,  the  model  at  hand  is  less  complex  than  the  one  dcvelcped 
by  Charnes,  Cooper,  and  Mellon.  Theirs  includes  a  maximum  vapor  pressure 
SDOclficatlon,  as  well  as  the  minimum  octane  condition  pp.  13fl-U;l  7. 

In  future  work,  the  concave  programminc  model  will  be  expanded  to  Incorporate 
additional  restrictions  of  bf  Is  nature.  Since  these  ether  properties  blend 
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llncarly,  they  can  b«  included  without  a  larpe  Increase  in  the  tbne 
required  for  computing. 

2i.  The  computational  procedure 

Kuhn  and  Tucker  hare  proved  the  following  "equivalence  theorem"! 

"Let  the  functions  fj^tx),  ...  f^(x),  g(x)  be  concave  ae 

well  as  differentiable  for  x  >  0.  Then  x°  is  a  solution 

of  the  maximum  problem  if,  and  only  if,  x®  and  some  u*^  give 

a  solution  of  t  e  saddle  value  problan  for  0(t,u)  S  r(x)  ♦ 

1/ 

u'Fx." 

For  the  gasoline  blending  case,  g(x)  is  given  by  equation  (3) 
above,  and  the  fj(x)  by  (6),  (7),  (9),  and  (10).  The  reader  can  verify 
for  himself  that  all  these  expressions  are  differentiable,  and  that  (3), 
(6),  and  (7)  are  concave.  Furthermore,  (9)  and  (10)  will  be  concave  over 
the  range  0  <  <  3.  In  otr.er  words,  if  we  have  the  answer  to  the 

problem  of  mlnlmaxlng  0(x,u),  we  also  have  the  answer  to  the  problem  of 
maximizing  (3),  subject  to  restraints  (L),  (5),  (6),  (7),  (9),  and  (10). 

y 

The  computational  problem  Is  that  of  finding  a  ainlmax  for  J5(x,u). 

1/  /~7,  p.  Ii86_7  The  expression  u'Fx  is  shorthand  for  the  following: 
uj^  fl(x)  ♦  U2  f2{x)  ♦...♦■  Uj^f^(x).  The  uj  are  non-negative 
Lagrangean  multipliers,  and  crrrespond  closely  to  the  economist's 
notion  of  "shadow  prices". 

?/  Note  that  in  the  present  version,  there  are  not  only  lower  limits  on 
~  the  Individual  x^  Imposed  by  expression  (L),  but  also  upper  limits  on 
the  X{^j ,  determined  by  (5).  For  computational  purposes,  It  was  easier 
to  impose  Uiese  upper  limits,  rather  than  to  set  up  additional  f«(x) 
restrictions.  The  upper  limits  make  no  essential  difference  in  the 
applicability  of  the  '^uhn-Tucktr  theorem. 
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The  computational  procedure  la  an  Iterative  one,  at  each  step  con¬ 
verting:  a  vector  x(t),  u(t)  into  a  new  vector  x(t  f  1)  and  uH  ♦  1).  The 
solutions  observed  have  tended  toward  the  minimax  point,  but  I  can  give  no 
strict  proof  of  the  necessity  of  this  convergence.'  The  new  vector  generated 
is  not  necessarily  "attainable",  nor  is  It  usually  "efficient"  in  the  sense 
of  Koopr.ans  £(  ,  p.  7?  7.  "Profite"  do  not  increase  monotonically  as  they 
do  at  each  step  of  Dantsig's  simplex  procedure  for  linear  progranmlng 
Despite  these  apparent  shortcomings,  the  algorithm  gives  useful  answers  in 
the  cases  examined  to  date. 

The  solution  must  be  started  off  from  some  Initial  vector  x(0),  u(0). 

In  principle,  this  may  be  any  arbitrary  non-negative  vector.  In  practice, 
though,  there  will  be  a  considerable  reduction  in  computing  time  If  a  "good" 
Initial  position  Is  selected.  At  each  step,  first  the  u  vector  Is  determined, 
an;  then  the  x  vector.  For  the  former,  the  basic  Iteration  consists  of  two 
steps,  and  for  the  latter,  of  four  steps.  The  problem  Is  one  of  determining 
A  uj  5  u(t  ♦  1)  -  u(t),  and  Ax^  S  Xj(t  ♦  1)  -  xi(t).  The  procedure  looks 
cumbersome,  but  less  than  two  minutes  of  machine  time  are  required  for  generat¬ 
ing  the  whole  set  of  A  uj  and  A  Xj^.  For  the  A  uj,  the  procedure  is  as  follows i 

1.  If  ^  ^  ^  0,  then  the  "candidate"  A  u.  =  -k* 

Sujft)  ”  ^  ^ 

/A 

If  ^  0  (t)  <  0,  then  the  "candidate"  A  u.  a  kj 

3uj(tJ 

Notei  kj  is  an  arbitrary  positive  constant. 


2.  If 

Uj(t)  ♦ 

Auj 

<  0,  then  A  uj 

II 

1 

c 

c* 

If 

Uj(t)  ♦ 

A  UJ 

^  0,  then  A  uj 

Z  Ujft) 

-  L,  .  -  -  -  -  w  -  ^  k  •  W  -  -  V  I.  w  a  J  .  w  .  W  w  W  .  a  .  la  W  W  W  tl  -'W  uX  _>  .V  .V 
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For  the  Ax^,  the  procedure  becoraest 


Ir  It,  for  any  j1,  fj(t)  <  0,  and  if  liLili  .  AElilLj  0»  thun 

(t)  dxi(t)  J 

‘  A  z  Oj  proceed  to  evaluate  A  ^  Otherwise,  proceed 
to  step  2. 

2.  If  ^  0  (t)  ^  0,  then  the  "candidate"  A  x^  r 

dxi  (t) 

If  d  Jf  (t)  <  0,  then  the  "candidate"  x*  r  -ki 

Note:  Is  an  arbitrary  positive  constant. 

3.  If  Xi't)  ♦  <  0,  then  A  a  -Xj^(t);  proceed  to  evaluate 

’H  e  1* 

If  Xi(t)  ♦  ^  ^  proceed  to  step  b. 

li.  If  ^  \  111,  and  if  1  ^  1,2,  A  -  A  x^;  proceed  to  evaluate 

*1  ♦  1- 

If  ^  =  Ll,  or  if  i  •  L2,  calculate  3  -  ^  *l.j* 

If  this  expression  is  non-nepative,  A  x^^j  r  A  X]^j .  Otherwise, 
A  Xj^  z  3  -  xj^j(t).  Proceed  next  to  evaluate  *i  *.  i* 

For  the  A  Uj,  th.e  explanation  of  this  ritual  is  straiptitforward. 

Step  (1)  tells  us  to  decrease  Uj  by  an  arbitrary  amount  kj  if,  for  a 
"s:;iall"  change  in  Uj  alone,  the  effect  will  be  to  decrease  0.  Similarly , 
there  is  to  be  an  increase  in  Uj  if,  for  e  ”sT.all"  Increase,  0  would 
decrease.  Step  (2)  prevents  Uj  from  becoming  negative. 
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For  the  A  the  Justification  la  more  tortuous.  (Steps  (3)  and  (L) 
are  the  obvious  ones  -  respectively,  lower  and  upper  limits  on  the  Indi¬ 
vidual  Xi.)  As  with  the  A  u«,  the  sign  of  S  0  (t)  is  primarily  the 

"3  *1  { t)  1/ 

.  criterion  that  determines  whether  to  take  a  ^losltlve  or  a  negative  step.” 
Preliminary  hand  calculations  Irvdlcated,  however,  that  the  unqualified 
'•direction  of  ascent  (descent)  rule"  would  lead  to  major  oscillations  during 
tne  process  of  convergence.  For  that  reason,  rule  (1)  ~  a  relaxation- 
type  principle  —  was  Inserted.  Translated  Into  English,  the  statenent 
reads,  "If  one  of  the  restraints  is  being  violated,  and  if  the  dlrection- 
of-ascent  rule  would  lead  to  an  even  greater  violation  of  this  restraint, 
then  the  particular  should  not  be  changed."  An  IBM  run  is  now  in  process 
for  purposes  of  comparing  the  efficiency  of  the  relaxation-type  algorithm 

y 

with  that  of  the  unmodified  di recti on-of-ascent  operation. 


1/  The  reader  should  note  that  0  is  being  minimized  with  respect  to  u,  and 
maximized  with  respect  to  x.  ” 

?/  Since  the  original  draft  of  this  manuscript,  the  run  has  been  completed. 

”  At  the  end  of  110  steps,  there  was  still  no  evidence  of  convergence. 

Parameter  values  and  initial  conditions  were  identical  with  those  for 
run  6.  The  only  difference  was  the  elLmlnation  of  step  (1)  for  the  computation 
of  A  xj^.  For  the  reader's  convenience,  the  time  series  on  three  of  the 
variables  aro  Indicited  below: 


t 

*31(1) 

Xp^Ct) 

Uj(t) 

0 

1,128.6 

571. Ij 

12.10 

10 

1,368.6 

651.1: 

11.90 

TO 

1,128.6 

151.1 

11.50 

30 

l,l.2H.6 

611.1 

11.90 

LO 

1,128.6 

11.90 

50 

1,308.6 

651.14 

11.50 

60 

1,165.6 

651. It 

11.90 

70 

1,2L6.6 

1491.1: 

11.30 

50 

1,‘<09.6 

65l.li 

12.30 

<30 

l,Oc8.6 

531. 14 

11.70 

100 

1,308.6 

731.6 

11.70 

]iU 

1,136  8. 6 

651. L 

11.70 

Table  1.  Fara-neters 


Run  Number 


St<  ps 
1-1L9 


Sbeps 

l50-3li^ 


r-, 


’  n 


;,i 


hr 


<  . 

y  *)> 


h.r 


y. 


k- 


12. OC 
0.23 
1 ,  JOu. 

io. 

10. 

0.1? 

LO. 

1.0. 

■^6. 

•J.l? 

0.(0 

O.io 

o.cu^ 

U.Ouf. 


ir.ou 

o.?3 

1,001. 


.  0 


r6. 


Ic. 

i;0. 

?6. 

0.  ( 
0.1‘^ 
J.l^' 

0 .0j3 
0.003 


ir.oo 

o.?3 

1,000. 

1.0. 

LO. 

0 


JO. 


1.0. 


L' 


0.1  •; 

U.1‘ 
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ir  Six  Sets  of  Calculations 


3 

li 

5 

6 

12.00 

ir.oo 

11.60 

12.00 

0.23 

0.23 

0.23 

U.60 

500. 

300. 

1 ,000. 

ro. 

ro. 

ro. 

<.  ^  • 

20. 

20. 

"0. 

('‘•ft 

30. 

■>0. 

30. 

30. 

o.O^ 

0.03 

0.03 

t'.03 

ro. 

20. 

?0  • 

f  '  • 

20. 

20. 

20. 

:j. 

"0. 

30. 

"0. 

•}  . 

V  a 

w.O"' 

0.03 

.  .03 

O.io 

0.10 

0.10 

0.10 

.'.10 

o.lo 

0 . 10 

j.in 

•  vV  V  ^ 

0  . 

.^O'.o 

0.-V3 

0 .  Ou  3 

o.O  '3 

.  .JO 3 

A*  the  preccnt  time,  I  am  in  an  ur.f  l<.asant  quandary  ccnc>.:nlnp  t..e 
loflc  of  these  conputatlonal  techniques.  On  the  one  hand,  I  r.avt  been 
unable  to  prove  the  necessity  of  convergence  for  cither  the  modified  or 
Uie  uinodlfiod  procedure.  On  the  other,  I  have  been  unable  to  construct 
a  cointer-cxample  for  a  case  where  the  and  the  kj  may  be  made  arbi¬ 
trarily  small.  Even  worse,  tl.ere  is  not  yet  any  general  criterion  for 
determining  whether  the  relaxaticn-type  procedure  Increases  or  decreases 
the  efficiency  In  ccmputlnp  out  a  particular  numerical  example.  Tnese 
are  all  problems  on  which  T  will  welcome  mathematical  advice, 

? .  Six  c-rputatlon  runs 

During  the  course  of  six  I*  rvms,  tne  only  nunerical  parameters 
altered  were;  pj^,  p-^,  the  kj^  and  the  kj.  The  individual  values 
for  these  are  listed  below  in  Table  1.  The  remainder  were  maintained 
at  the  values  indicated  in  Appendix  I. 

In  order  to  provide  the  most  severe  test  of  the  computational  pro¬ 
cedure,  the  first  run  was  made  with  a  distinctly  non-optlmal  initial 
position  —  in  fact,  from  the  zero  vector  for  all  xj.  (Exact  numerical 
values  for  Uio  x(0)  and  u(0)  of  each  run  are  listed  in  Appendix  II.) 

This  non-optlmal  feature  shows  up  in  the  time  series  for  the  first  20 
steps  (Kigurer  3-7).  Most  of  the  x^^  Indicate  a  strong  upward  trend 
during  this  Initial  phase.  The  only  legitimate  inference  tnat  can  be 
drawn  from  the  series  at  step  (20)  is  that  nothing  like  a  stationary 
equilibrium  has  been  readied.  In  fact,  there  continued  to  be  signifi¬ 
cant  trends  in  the  x^  until  r0(j-300  iterations  naj  been  performed.  The 
process  wus  linally  stojped  at  the  3li6tii  step.  The  series  for  steps 
arc  reproduced  as  Figures  3-12,  and  the  vector  at  this  point 
Is  given  in  Column  A  of  Table  2, 


V'? 
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Even  at  the  3li6th  step,  the  systfsi  Is  not  at  a  5ad<!le  point.-  According 
to  the  graphical  solution  (column  C),  at  the  optlnrum  i>oint  r  xj2  = 

^22  ^  *32  ”  ^i^ough  the  IBU  tabulation  does  indicate  zero  levels  for 
X2^2  *32»  nevertheless,  Xj^j^  and  x^j  both  show  up  at  low  positive  levels. 

More  serious,  the  exact  solution  is  infeasible.  For  example,  the  schedule 
requires  more  of  stock  2  than  is  available  —  1,026.1  gallons,  as  against 
the  1,000  galloit  limit.  Cnl/  a  simple  adjustment  is  required,  though,  in 
order  to  obtain  feasibility.  ''The  entries  for  Xjj^,  X31,  and  xjj  were 
each  multiplied  by  1,000.0  s  O.bTLfS.  In  order  to  satisfy  ft,  it  was 
necessary  to  make  a  slight  increase  In  x|ji{  f|^  could  be  satisfied  with  a 
lower  value  of  finally,  it  was  noted  that  the  adjuctnent  to  the  x^  left 

li3.6  gallons  of  stock  1  in  excess.  Since  stock  1  can  be  sold  as  premium 
grade  gasoline  by  adding  only  0.60  ml.  of  lend,  this  li3.6  gallon  excess  was 
credited  to  profits  a£  a  value  of  ll«dl6  ^/gallon.  The  adjusted  blending 
schedule  is  both  feasible  and  efficient.’  The  payoff,  g(x),  amounts  to  $230.53 
as  against  a  possible  S231.03,  indicated  by  the  graphical  solution.  These 
two  nucbers  coincide  to  within  1  part  in  500  —  well  within  the  limits  of 
accuracy  of  the  graphical  method. 

Having  established  the  general  convergent  nature  of  the  alporitfr:,  the 
remaining  calculations  were  carried  out  froa  Initial  vectors  that  were  con¬ 
sidered  to  be  sensible  starting  point? .  Run  2  was  started  off  fron  C.,  the 
vector  that  is  labelled  "graphical  solution"  in  Table  2.  The  solution  never 
departed  markedly  from  Uie  initial  values,  and  was  halted  at  the  63I^1  step, 
giving  a  second  column  labelled  "A.,  raw  vector"  in  Table  2.  The  rough 
solution  was  adjusted  in  the  same  way  as  has  been  described  for  run  1,  and 
the  results  also  entered  in  Table  2.  '.'he  adjusted  g(x)  --  ?2?3.62  —  is 
araln  acceptably  close  to  that  for  the  graphical  solution,  *223,73. 
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10 

Step  number 


Fig.  5— Steps  1-20  of  run  number  I 
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Fig.  B  — Final  20  steps  of  run  number  I 
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Fig.  9  — FInol  20  steps  of  run  number  I 


Fig.  10— Final  20  steps  of  run  number  I 
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Table  2.  Solution  Vectors,  Runs  1  and  2 


Run  1 


Run  2 


A. 

Htw 

Vector 


B. 

Adjusted 

Vector 


Craphic  U  Ra« 


B. 

Adjust»*:l 


Solution  H  Vector  I  Vector 


C. 

''inphlc 

Holution 


final  L 


1,285 


5L5.6 

l,251i.'> 


1,L25.6 


] 

l,L)OC'. 


fl(») 

f-fy) 

e(x)  (  j  illara)  231 .SI 


'X).55  231.03 


3  J.HL 


2?3.<55  2?8/>2  ?23.73 


383 

ri- 


Table  3*  Scl'itlon  Vectors, 


rt'in  3 

A 

Raw 

Vteter 

B 

Adjusted 

Vector 

G 

Cri-phic 

Final  t 

70 

~ 

— 

*1] 

1 

0 

0 

1 

! 

*?1 

705.7 

601. q 

nt 

>^31 

71J  .3 

'/CK\3 

ni 

"'ll 

2.79 

2.f0 

2.75 

0 

0 

0 

*r? 

100 

qd.o 

70 

CM 

K 

0 

0 

0 

0.75: 

1 

0.7^ 

0.7U 

fl(x) 

1 

o.e 

1 

fpCx) 

— 

0 

— 

>,(x)  ;  i  ] 

ari.)  172.1? 

Kq.'JB 

170. 11 

A 

Vecter 

30 

C 

r3i.i 

L2-.6 

3.0 

0 

■<80 

0 

0.79 

I 


A'lju-.ti  .1 
V'fcbci- 


0 


uV  A. 


3.0 

0 

I4OO 

0 

0.70 

0 

0 


110.69 


IIS.15 
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Table  li.  Solution  Vectors,  Runs  5  and  6 


Run  5 

[ 

Run  6 

A 

B 

C 

1  ^ 

B 

C 

Raw 

Adjusted 

Raw 

Adjuote  1 

Orai.nic 

Graphic 

Vector 

Voctor 

Solution 

Vector 

Vector 

Solution 

Final  t 

80 

— 

eo 

— 

— 

"ll 

21i0 

233.0 

2^7 

80 

70.6 

156 

’'a 

LO 

38.5 

0 

0 

0 

0 

^31 

1,068.6 

1,037.5 

1,090 

1,158.6 

1,1?9.2 

i,l31i 

1.05 

1.15 

1.05 

1.02 

1.05 

1.10 

*12 

0 

0 

0 

0 

0 

0 

*2? 

591. li 

571.. li 

673 

611. <4 

601.? 

6,98.3 

6o 

58.2 

0 

180 

177.0 

21.7 

*12 

0.55 

0.58 

0.70 

0.1i6 

0.50 

0.66 

— 

0 

— 

— 

0 

— 

f^^x) 

— 

58.3 

— 

— 

L.O 

— 

»'(x)  (dol 

1 

larc  )r2l,  .01 

223. L5 

225.9? 

222.15 

2.1.91; 

V) 
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Durlng  runs  1  and  2,  the  rraphical  solutions  mrc  lettrmlned  In 
advance.  For  3-6,  however,  they  were  obtained  after  the  adjusted 
solution  had  been  calculated.  (For  runs  3  and  h,  the  initial  ;:utss 
turned  nut  to  be  the*  solution  to  the  problem.)  Table?  3  and  L  reveal 
tnat  the  algorithm  continued  to  perform  In  a  satisfactory  way  during 
the  last  four  runs.  For  run  6,  in  fact,  the  adjusted  solution  at  step 
'30  inilcates  a  higher  net  realisation  than  the  graphical  technique  — 
Crft.l^,  as  against  ^PFl.'Jli. 

6.  Interpretation  of  results 

T:ie  nunerlcal  parameters  used  in  run  1  for  P2^,  p->,  and  qo  were  con¬ 
sidered  to  be  the  basic  set,  and  In  each  subsequent  run  not  more  than 
one  of  these  was  altered.  In  (♦'3  and  ^It,  the  value  of  q^was  changed;  In 
run  5»  P^i  and  in  run  6,  p^.  The  adjusted  results  for  all  the  six  runs 
are  summarized  by  Table  5  In  terns  of  four  indicators  —  the  output  of 
premium  crade  gasoline,  the  output  of  regular  grade,  the  input  of  TtL, 
and  the  refinery  realization.  These  data  are  also  presented  on  Charts 
13-16. 

Run  2  differed  frooifj'lin  only  one  significant  respect.  IHirin^;  it, 
tne  refiner  was  not  permitted  to  make  stock  3  —  the  fictitious  inter¬ 
mediate  proluct.  The  initial  olending  materials  had  to  be  assigned 
directly  to  the  finished  products.  In  other  words,  the  refiner  was  not 
penniLlt-d  to  take  full  advantage  of  the  concave  octane  nu.aber  blending 
relationship. 

This  additional  restriction  had  no  effect  upon  the  general  strategy 
of  converting  virtuaDy  the  entire  volume  of  blending  materials  into 
premium  grade  gasoline.  The  restriction  ^e«  have  the  effect  of  increas¬ 
ing  the  TEL  requirements,  and  of  decreeing  the  refiner's  not  realizatinn. 
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"  ht  dlfferenc*'  In  realization  at  first  gl-ance  appears  snail  —  S22‘’.6r,  as 

iiialnst  •?30.5tl  in  run  1.  However,  for  a  refinery  turning  out  20,000 

barrels  of  gasoline  per  calendar  day,  such  a  discrepancy  would  come  to  an 

annual  figure  of  roughly  300, 000*.  To  put  the  matter  another  way,  the 

flrat  blending  schedule  does  the  same  Job  as  tl.e  second,  with  a  total  Trl 

r*-ip>irer.ent  that  is  about  20?  smaller.  From  this,  I  would  conclude  that 

It  Is  eminently  worthwhile  for  a  refiner  to  incur  consi itrable  expense  in 

testing  gasoline  blends,  and  In  attempting  to  predict  the  occurrence  of 

curvilinear  octane  relationships. 

In  the  first  run,  it  turned  out  tint  the  most  profitable  blending 

schedule  was  ono  that  assigned  all  2,-00  gallons  Into  premi-um  grade  gasoline. 

This  result  was  quite  in  line  with  a  rule-of-thinb  that  is  c  xurnonly  heard 

among  refiners  --  to  wit,  that  the  max'mum  profits  will  be  obtained  by  making 

1/ 

the  oiaximum  amount  of  the  high-priced  product.  This  rule-of-thunb  is  not 
shaken  by  the  results  of  run  3  or  h.  At  a  availability  level  of  500 
gallons,  the  best  strategy  is  to  produce  virtually  nothing  but  premlun  grade 
gasoline.  At  a  qj^  Itvel  of  ‘’00  gallons,  however,  tne  3t0  ml,  limit  on  TFX 
concentration  has  cone  into  effect.  A*,  this  lew  level,  the  best  the  re¬ 
finer  can  do  is  to  add  3.0  ml.  to  each  gallon  of  priniuu  grade  gasoline,  and 
and  to  assign  the  left-over  stock  2  to  ’he  proiictlon  of  the  regular  grade 
product . 

The  preceding  results  dc  nat  speak  wi^ll  for  an  alternstlvoly-hcard  dictum 
that  the  best  general  course  to  follow  is  to  assign  3.0  ml.  cf  ethyl  fluid 

1/  I  do  not  allege  tfiat  refirars  generally  follow  such  a  rule  —  only  that 
”  they  say  LJ  oy  do. 
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9 

a 

c 


Fig.  13 


Net  realization  (dollors) 


to  cic.*’!  1'alic.r.  'f  prolact.  Only  in  -naVlng  prfmium  graJe  pasolir/.  during 
r’an  I  did  It  fnv  to  use  t^e  m  ixlunax  allowable  TEl  concent  ration. 

iluns  ^  and  6  cast  doubt,  not  only  on  the  3*0  ml.  principle,  but  also 
on  the  T.axlm'iE-pro Juctlon-cf-prc'Tilu.n  theory.  In  run  5,  the  only  para- 
Tiettr  ci'.ange  from  *1  was  e  reduction  in  the  price  of  premium  grade 
pr-oduct  from  1  ?,*  'ja'.  1  on  to  11.0.  T-i-  difference  in  tho  price  structure 
is  sufficient  to  bring  about  a  substantial  increase  in  the  optimal  out¬ 
put  of  regular  grade  product  —  from  less  than  130  gallons  in  ffl,  to 
aLr.ost  700  gallons  in  Tr.ir  shift  is  also  accompanloi  by  a  marked 

reduction  in  the  total  Input  of  ethyl  fluid. 

For  run  6,  a  different  sort  of  change  in  the  urice  structure  was 
tried.  T:  e  pr^.ce  of  ethyl  fluid  wci  virtually  trebled  —  from  0.r3d/nl. 
to  0.{'0i.  Again  comparing  tne  results  with  there  is  a  shift  into 
the  production  cl  regular  graJo  j-asollne,  and  a  reluction  in  tt.e  input 
of  c»'.y]  fluid. 

To  put  the  matter  another  way,  the  rulcs-of-thimb  imply  tnat  the 
same  input  and  output  rt  1  ■.  11  onr.hip  s  are  op  timal  un  ler  all  price  struc¬ 
tures.  One  of  these  rules  does  iiappen  to  t  iv  a  good  pc  rfonr.ance  for 
a  p.irtlcular  c.rrint  set  of  r:-ic*n,  hut  1  *^6  not  apply  outside  this 
li-'tted  rui’e,  ''nfort’ir.ately ,  tb.f're  U'  no  <  ijy  rseape.  li.e'  only  way 
to  v<-rify  one  of  ♦hes’--  r ul  “r: -of-thumb  in  a  concrete  ca'-e  is  to  carry 
th^ou a  Ic’ailo]  I'or'oal  analysis. 


/.  I  u  cions 

At  this  foinl,  it  is  a:  ;  roprlatc  to  cum  up  the  findings.  Or.  ♦  h'O 
n'gfiv"  i:  li  ,  the  study  h-,s  demonstrated  that  certain  rule  s-of-thamb 
worV  : -iJly  —  even  ♦>hongh  these  principles  have  occasionally  been  put 
lorwiri  hy  -ipn  in  t/iO  refining  iniurtry.  The  s^.udy,  moreover,  nas  mown 


.  -v*  f\M  rje  "ji 
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that  the  oftimal  input  and  output  relationships  will  be  intLnately  linked 
with  the  price  structure.  Figures  13  and  llj  give  some  clues  as  to  the 
nature  of  thito  linkage.  Before  the  economist  jumps  to  the  conclusion  that 
these  really  arc  supply  and  demand  curves,  however,  he  is  urged  to  reread 
the  caveats  of  Section  3. 

Cn  the  positive  side,  the  claims  are  necessarily  more  cautious.  The 
study  indicates  that  a  refiner  will  do  well  to  make  a  careful  inver.ti^'ation 
of  his  gasoline  stocks  for  non-linear  blending  characteristics.  The  results 
do  not  imply  that  a  20f  saving  in  TEL  costs  will  inevitably  occur. 

Wore  significant  from  a  long-range  viewpoint,  the  computing  runs  indicate 
that  it  is  possible  to  solve  concave  programing  problems  on  a  digital  com¬ 
puter.  The  particular  procedure  employed  is  undoubtedly  wasteful  of  machine 
time,  and  superior  methods  will  surely  evolve,  Tne  next  programming  problem 
tnat  will  be  attacked  is  one  that  has  been  formulated  by  one  of  the  major 
dest  Coast  refiners.  Although  this  coming  problem  has  ?2  and  8  Uj,  the 

basic  notion  is  the  same.  The  quantities  available  of  various  blending 
stocks  will  be  taken  for  granted,  and  the  analysis  carried  on  from  that 
point.  Certain  cutting  temperatures  will  also  enter  as  Independent  vari¬ 
ables.  At  a  later  date,  it  is  expected  that  the  model  will  be  expanded  to 
incorporate  such  process  variables  as  reactor  temperatures  and  recycle  ratios. 
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Appeixilx  I.  Constants  f^r  Sly  Sets  of  Calculations 


P-3^3 

- 


p,  r  12^/,-il.  ) 

) 

) 

Ni  -  "^O  octane  ) 

number  ) 

) 


Sources: 

Ref,  8  i^lvcs  the  follcwlnp  Gulf  Coasti 
bulk  carcocs  price  range  for  90  octant 
(ASTM  Research)  premium  motor  gaeolinet 
n.7s:  -  1?  -  IP.r^tf/gal.  These  prices 
were  In  effect  cn  November  PL, 


P2  =  ll//ral.  ) 

) 

N’t  :  8"^  octane  ) 

number  ) 

) 


Hef.  8  gives  t.he  fcllcwing  Gulf  Coast, 
bulk  cargoes  price  range  for  t}  '^clar.e 
(kSTV.  Hesearch)  regular  motor  gasoline: 
10.7c  -  11  -  ll,?5vVgal.  Thoct  r-rices 
wert'  'n  effect  on  Novcr,ber  Oi, , 


I''  z  of  ) 

Til  ) 

) 


The  Loo  Angeles  Silos  Cffiei  of  *.hf  Lthyl 
Cerf'Orition  gave  ♦hie  quotaticn  rn  Nevemher 
Cl,  1^52,  It  applies  tc  motor  mix  eb.^vl 
fluid,  delivcrt'l  in  tank-car  l"'ts  to  my 
U,  b.  jestl  tialicn. 


r  ■<00,  ^00,  1,000  gal.  ) 

Ref.  9, 

of  catalylic-cr- ckcJ- ) 

irwn  of 

clr.s-golymcr  gasoline) 

> 

.’enuary 

;  i,0^)j  gol.  of  td.criial) 

the’T.sl 

crocked  r.i.'':line  ) 

cat  a.y’ , 

) 

day. 

p#  inJi'  .tcs  too  fclicwlnt' 

U,  £.  ;r'icl-.in,j  {.lant  csjacit’'  of 

1,  l<3f?: 

cracked:  >  arrvlr/  iay i 

rr.icke'l  jasalin*:  t:a:it?s/ 


.‘1,3.  There  arc  wicie  vuriaticrii.  in  t;.r 
ratio  uf  thermal  to  cata.’ytic  capicity 
as  between  in  iiviJual  rcfi  ier^ts  an.i 
refining  r»^gio.is. 


A{jpendix  T.  (continued) 


The  hj,  b,,  aM  Ci  coelficlents  for  eeel.  stock  .er.  derived  froci 
t„e  boeen  .nd  Nichols  d.t,  ^1-  P- 

eels  of  observations  are  shown  on  figure  1.  The  Individual  coef 


ficients  had  tue  foilcvtin^  numerical  values: 


Stock  n 


d'5.6 


-ir,]; 


9h.l 
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Appendix  II.  Initial  Vtctorc  for  Six  Sets  of  Calculation 


X)  J/ 

x^^(O) 

w  t 

^32'°) 

X|j2(0) 
Uj^O) 
UpIU) 
u-jf  J) 

Uj^(G) 


1,1)00.0 

l.CiOO.O 


?.?0 


i.'.oo  ir.io 

1C.  it'  io.7!r 

U.lOO  O.l'jK 

0.1'lJ  o,'.  \.l 


0.7 

.10  ir.o 


5  ! 

6 

1 

0 

0 

v) 

1,1:2?  .6 

•  1  /, 

^  1  ,  L.  t  .  . 

1.? 

1.2 

0 

0 

571.1: 

‘■71.1; 

0 

0 

0.7 

0.7 

12.10 

12. 1. 

10.75 

10.75 

0.1 5I4 

0.151; 

0.15L 

0.15-li 

Ajt1d«nduin  to  P-383 1  "ConcaTe  TYogranning  for  Oiiollrw?  Rlorulo 


In  c)l8cu*«lng  r-383.  Dr.  >jT.rtln  boctaian  uucoverad  a  problem  to  vblch 
I  had  not  ptavlouBly  given  attention.  He  haa  made  the  p-olnt  that  for  n 
function  of  many  variable*,  f(x),  to  be  truly  concave,  It  1*  necesoory 
but  not  Buffi  dent  for  that  function  to  be  concave  In  each  variable  taken 

■^parately . 

The  baalc  problem  arlae*  out  of  the  fixet  that  the  function  xy  la 
neither  concave  nor  convex  -  for  nonnagatlve  Independent  varlableo  x  and 
y.  Thi*  la  de*plte  the  fact  that  It  la  concave  and  convex  In  either  variable 
separately. 

Proof:  Conalder  tvo  point*  (x^^,  and  (r^,  y^)  Convexity  would 
require ; 

L  ’>•’>  ]  ['>•'>]  = 

^  ^  *2  ^1  -  *1  ^1  ^  *2  ^2 

(2)  x^  (y^  -  Xp  (y^  -  y^) 

For  nonnegntlre  x  and  y,  and  for  y^  ^  expreanlcn  (2)  requires  x^ 

<  x^.  But  the  cixivexlty  condition  ajTplie*  to  any  arbitrary  (x^.  y^)  aivi 

(xg,  3^2^’  *  oontradlctlou.  Slallarly  for  concavity.  Q.E.D. 

Rov,  Vhat  art  the  Ijoplloatlona  of  thlo  retult  for  the  gaaolliM  blending 

rolel?  F^rtumtely,  It  •  irne  out  that  by  a  re-deflnlllon  of  variables  .  the  pay¬ 
off  function  g(x)  can  be  made  concave  --  Indeed  llneor.  Unfortunately .  thimgh . 
thl*  re-def Inlt  1  on  d  leB  not  lead  to  any  rlgorou*  proof  of  the  concavity  of  reotmtnt 
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eqjallr^nB  (9)  ajid  (lO).  the  contrary,  It  shovs  up  the  pcsBlblllty  of 
non-conrA/lty  In  these  restraint  coodltlons  when  two  of  the  bleivllng  stocks 
have  widely  divergent  lead  Busceptibllltles . 

Let  ua  define  new  variables,  x.^  and  In  the  following  way: 


•  ’'21  *  *31’ 

’'32  **12^2  *■  *22  *  *32* 

Equation  (3)  can  nc?w  be  transfcmiedi  Into  a  linear  fom  involving  the 
tw  new  variables,  x^^  and  x^^  * 

(3a)  g(x)  Pi  .  Xj^  *  Xj^)  ‘  Pj  (X;lj  •  Xpj  *  *32*  '  P3(*5i  *  *52) 

To  simplify  the  discussion  of  equations  (9)  and  (lO),  I  aholl  arbitrarily 
assuiT^  that  both  x^^^  and  arc  constrained  to  be  terc,  and  will  focus 
attention  on  the  function  f^  (x),  dealing  with  the  octane  specification  of 
the  p.-emiuffl  grade  gasoline.  Suppose  that  the  octane  nun4>er8  and  the  quantities 
available  of  the  two  blending  stocks  had  been  at  follows: 


Stock  1 

stock  ;? 

Cuontlty  available, 

1,000  gal. 

l.OCO  gal 

Orlnne  niiraber,  clear 

■!*.0 

<X).0 

Octoine  number,  with  3  tfd. .  cf  TEL 

<0.0 

Namber  of  ml. /gal.  required  for 

X)  octane  leaded  gasoline, 

3.0 

• 

0 

Total  number  of  ml.  required  for 
X)  octanr  lrade*1  gntollne  ,  x,.^^ 

3,000 

0 

Assuming  straight  line  Interpolation  at  constant  lend  level.*)  for  the 
two  stocks .  and  assuming  a  SO-SO  blend,  we  come  up  with  these  resulto  for 


the  blend; 
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Q'jantlty  BLVullable,  (x^  +  *21^  2,(XXD 

Oc-t«jie  number,  clear  82.0 

Octane  number,  with  3  »!•  of  TEL  92.h 

Number  of  ml. /gal.  required  for 
90  octane  leaded  gasoline,  1.6 

Total  number  of  ml.  required  for 
90  octane  leaded  gasoline,  x_.  3*200 


In  other  worda.  If  the  two  initial  materlala  are  ethyl l»ei  «mJ  marketed 
separately,  the  average  TEL  requiremant  per  gallon  comes  to  1.5  irL  .  If  the 
two  are  blended  together,  the  same  amount  of  premium  gasoline  la  produced, 
but  the  average  lead  requirement  increases  to  1.6  ml.  Tills  behavior  clearly 
violates  concavity  of  the  function  f^  (x^ ^ *51^'  lends  to  the 
pooalhlllty  of  multiple  Isolated  nsxiiaa.  Lacking  the  eoncavlty  property  in 
our  functloni',  we  may  be  led  to  adopt  blending  schedules  that  are  locally 
optimal,  but  which  do  not.  In  fact,  represent  the  best  of  all  solutions  avail¬ 
able  . 

Nov,  how  serious  are  the  possibilities  of  such  local  aaxlnn?  The 
problem  cannot  be  shrugged  off,  but  I  believe  that  in  most  particular 

f 

applications  these  local  extrema  can  be  detected  by  th<;  Investigator.  In 
order  for  the  violation  of  concavity  to  take  place,  soar  of  tlie  blending 
material  must  be  particularly  high  la  lead  susceptibility,  and  sane  of  It 
must  be  extremely  lew.  In  addition,  the  first  imterlal  ^st  require  a  high 
TEL  concentration  In  order  to  meet  the  specification,  and  the  second  n  low 
TETi  level  Enrvlng  something  shout  the  occurrence  of  such  pitfalls,  the 
Intelligent  refiner  sbcruld  be  «ible  to  avoid  ensnaremnnt . 

The  airrieultles  arising  out  of  non-con^ uvlty  are  not  unique  to  the 
programming  method  described  in  r-383  The  sane  argument  would  aprply  if 
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and  vero  taken  as  parameters,  amd  a  separate  linear  programnlng 
problea  then  non  off  for  each  of  many  conblnatlons  of  the  two.  This  dif- 
flcTilty  would  arise  In  the  case  of  any  triad. -and -error  method  that  consisted 
In  changing  ?ne  variable  nt  a  time,  and  observing  the  effects  upon  the  payoff 
g(x).  And  finally,  such  failure  of  the  proper  curvature  conditions  could 
frustrate  an  optimliatlon  via  a  Lange-Lemer  shadow  price  market  nechanlsm. 

The  dllerma  raised  by  Dr.  Beckman  ahoxild  not  be  viewed  as  an  absolute 
Impasse,  but  rather  as  a  challenge  to  the  Ingenuity  of  workers  In  this  field. 
These  latter  can  take  comfort  In  recalling  that  truly  rigorous  analysis  Is 
seldom  applicable  within  the  domain  of  applied  npAhemtlcs. 


K.'y  V.  V 
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Mr.  Lloyd  Phlllpaon  haa  pointed  out  the  follcviixg  misprints: 


last  line  of  p.  12  should  read: 

line  6  of  p.  13  should  read:  ** 
/\ 


"If  u  (t)  +  A  u  >  0,  then  A  u  ^  A  u.(t)" 

0  V  u  J 


<  0,  then  the  'caivl  I  late* 

axj  u) 

r 


w  U-*  U-H  L'>  L*>(  1/^  jTX  ■ 


